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Abstract 

o 

The two-step model with tt, p and u exchanges taken into account is applied to 
investigate the reactions pn — > dr] and pn — > dn' near the corresponding thresholds. 
The existing experimental data on the reaction pn — > dn are analyzed and predic- 
tions for the cross section of the reaction pn — > <i?/ are presented. It is found that 
7T- as well as /9-meson exchanges yield significant contributions in both reactions. 
Furthermore, the effect of the rjN final-state interaction is studied. It is shown that 
an r/iV-scattering length of Re a v ]y < 0.3 fm is sufficient to reproduce the observed 
^ ■ enhancement of the pn — > dn cross section close to threshold. 

PACS 25.10.+s; 13.75.-n 
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Recently, the total cross section of the reaction pn — > dr] has been measured 
close to threshold at CELSIUS in Uppsala [1,2]. These data are remarkable for 
a variety of reasons: E. g., they reveal that an older estimate yielding a rather 
large near-threshold production cross section for the reaction pn —>■ dr) [3,4] 
is definitely ruled out [2]. Furthermore, they show that the production cross 
section very close to threshold is significantly enhanced over the expectation 
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for a two-body phase-space behaviour [2], indicating the presence of a strong 
final-state interaction. 

Though there is a large body of theoretical papers about the reaction pp — > 
ppi] only a few of them have also dealt with aspects of the near-threshold 
pn — > dn and pn (pn)i=on production [3-7]. Moreover, all of them rely 
on the old and now obsolete data point of Ref. [3]. Thus, in the following 
we want to present a timely re-analysis of this reaction utilizing now the 
new data from CELSIUS. In particular, we are going to present results of a 
microscopic calculation of the pn — > dn cross section within the framework 
of the two-step model (TSM), previously applied to the description of the 
Pontecorvo reactions pd — > pM (see, e.g., Ref. [8] and references therein). We 
start from the n exchange contribution to the 77-meson production and study 
the relevance of other production mechanisms such as the p and uj exchange. 
Furthermore, we investigate the effect of the final-state interaction between 
the r) meson and the nucleons employing an optical-potential model as well 
as a multiple scattering approach. Finally, as an application of our model, 
we present predictions for the cross section of the reaction pn — > dn' near 
threshold. 

In the TSM the reaction pn — > dM proceeds in two steps: i) the initial nucle- 
ons produce an 77 or 7/ meson (denoted as M) via virtual ir-, p- or w-meson 
exchange; ii) the proton and neutron form a deuteron in the final state (see 
Fig. 1). 

The amplitude which describes the contribution of the 7r° exchange can be 
written in the form 

Tpn^dui 8 ! t, u) = Ap n _^ dM (s, t) + Ap n ^ dM (s, u) , (1) 

where M is the pseudoscalar meson 77 or n' . s = (pi + P2) 2 , t = (P3 — Pi) 2 , 
and u = (p3 — P2) 2 where pi, P2, P3, and P4 are the 4- momenta of the proton, 
neutron, M-meson and deuteron, respectively. Treating the nucleons inside the 
deuteron nonrelativistically, the two terms on the right hand side of Eq. (1) 
are given by (see, e.g., Ref. [9]) 
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where e d is the polarization vector of the deuteron and tp\ the spinors of the 
nucleons in the initial state. f w is the nNN coupling constant, for which we 
use the value f%/(4ir) = 0.08. Furthermore we have 



M>1) =V2^J (£ + pi) -^j- 2 , (4) 

**(*,P1)=-T^T» (5) 
with the kinematics 



q 2 = ml - 5 (k 2 + /3(pi)) - IpS, q = k + pi , 
P(pi)=(p? + rnl-T?)/6 , S = 1 + 71/m, ^ = yfft + m 2 - m . 



In Eqs. (4) and (5) ^(/c) is the deuteron wave function, F n (q 2 ) the form factor 
at the 7rNN vertex and m the nucleon mass. Note that the amplitude M n 
corresponds to the exchange of a neutral n meson only (see the left diagrams 
in Fig. 1). To take into account also the charged pion exchange we have to 
multiply amplitude (1) with an isospin factor (which is equal to 3). 

The form factor F w (q 2 ) is taken to be of monopole type with a cutoff mass 
A,r = 0.8 GeV/c, in line with recent QCD lattice calculations [10] and other 
information [11]. The D-wave part of the deuteron wave function gives only 
a small contribution and can be savely neglected. This was demonstrated 
in Ref. [8] for the case of the reaction pd — > Mn where the same struc- 
ture integrals (4) for 7r, p and uj exchanges occur. Note that the property 

— * — * 

M 7r (— pi) = — M 7r (pi) ensures the correct behaviour of the amplitude (1) un- 
der permutations of the initial nucleons. The amplitude should be symmetric 
in the isoscalar state. 

Near threshold only the S-wave part of the amplitude of the elementary reac- 
tion irN — > MN should be of relevance for the production amplitude T pn ^ dM . 
It can be obtained from the experimental S-wave cross section through the 
relation 

pcm 

\A 7T o N ^ MN (s 1 ,t)\ — \A n o N ^ MN (s 1 ,u)\ = 8ns 1 ^^a n - p ^ Mn (6) 

Pm 



where S\ is the square of the cm. energy in the Mn system. The values of 
these S-wave cross sections that we used in our calculation are 
<7«-„ = (21.2 ± 1.8)p c ^ fib and ^- p ^ n = (0.35 ± 0.03)p c ^ /ib 
(p^ in MeV/c). They are extracted from the data given in Ref. [12]. 
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The total cross sections resulting from the 7r-exchange contribution (with = 
0.8 GeV/c) are shown in Figs. 2 and 3 (dashed curves) as a function of the 
cm. excess energy Q = y/s — m 3 — m 4 . The data points for r) production 
(Fig. 2) are taken from [1,2]. 

Evidently pion exchange alone cannot describe the 77-production cross-section 
data. Of course, we could have simply adjusted the cutoff mass in the pion 
form factor in such a way that we reproduce the data. But it goes without 
saying that the exchange of heavier mesons should also be considered as well 
as effects of the initial- and final-state interactions. 

In order to obtain a direct estimate of the contributions from the p and uj 
exchanges we calculated them in the framework of the TSM using the vector- 
dominance model (VDM) prediction for the amplitude pN — > rj(rj')N and 
assuming that Aun-^^n ~ A p o N ^, v r v /\ N . We derive the S-wave 7 AT — > r)N 
amplitude from the Mainz data at E 1 = 716 MeV [13] where the differential 
cross section of the reaction ^p — > rjp is isotropic. In the case of rf we take the 
7P — > rj'p cross section of about 1 - 1.5 pb at E 1 = 1.5 GeV from Ref. [14]. 
According to Ref. [14] the differential cross section of the reaction 7p — > rj'p at 
this energy is almost isotropic. In the actual evaluation of the amplitudes we 
use the coupling constants 7 P and 7^ of Ref. [15]. The VN — > MN amplitudes 
for longitudinal polarization of the vector mesons V are taken to be the same 
as the transversal ones. Defining the cross section averaged over the V-meson 
polarizations as dcr/dfi = (p/ m /pf m )|/| 2 we obtain the values \f\ 2 = 720 /xb/sr 
and 60 - 90 pb/sr for the reactions p°p — > rjp and p°p — > r/'p, respectively. 
Those numbers should be compared with 360 and 9.8 pb/sr for the reactions 
7r°p — > rjp and n°p — ■> rj'p, respectively. 

The amplitudes for the vector-meson exchanges can be written in the form 
(treating again the nucleons in the deuteron nonrelativistically) 



G 

A ^dM(s,t) = ^(pl 2 (p 2 )(-i cr 2 ) • A v o N ^ MN (s 1: t) x 
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and 



M^(p 1 ) = V2^ f[(k-p 1 ) + ^(k+fc)} x 

J /77 t / 




GV is the vector coupling constant and «y is the ratio of the tensor and 
vector couplings. For these coupling constants and for the corresponding vertex 
form factors we use the values from the full Bonn AW potential [16], i.e. 
G 2 p /4n = 0.84, k p = 6.1, G 2 /4vr = 20, = and A p = 1.4 GeV/c, A w = 
1.5 GeV/c. Note that the amplitudes (7)-(8) should also be multiplied by the 
isospin factor 3 in the case of p exchange. 

As far as the initial state interaction (ISI) is concerned it would be, of course, 
desirable to calculate its effect by means of a realistic model of the AW interac- 
tion. However, there are basically no AW potentials available in the literature 
that are still valid in the energy range where 77 and, in particular, 7/ production 
occurs. Therefore we decided to introduce an overall normalization constant 
which effectively accounts for the expected reduction of the production cross 
sections due to the ISI. This constant is adjusted to the data on the reaction 
pn — > dn and then used for calculating predictions for the reaction pn — ► dn'. 
The latter step is based on the reasonable assumption that the ISI does not 
change that much in the energy range from the n- to the ^'-production thresh- 



Alternatively, we employ an estimate for the effect of the ISI that has been 
presented recently in Ref. [17]. It is based on the assumption that, at the large 
energies required for the production of heavier mesons, the AW t-matrices are 
basically constant for (off-shell) momenta not too far away from the on-shell 
point. In this case the effect of the ISI can be expressed by a simple reduction 
factor A that depends only on the AW phase shifts and inelasticities, i.e. [17] 



In our case the initial pn system is in the 1 P\ state. The corresponding phase 
shifts and inelasticity parameters near the thresholds of n and n' production 
can be taken from the VA phase-shift analysis [18]. The values at 7i a b ~ 1250 
(2400) MeV are 5 L (pi) = -29°(-2.1°) and n L {p{) = 0.54(0.828). This yields 
a reduction factor of Aisi = 0.46(0.83). 



old. 



KTnn) = Vl(T nn ) cos 2 (6 l (T nn )) + 1/4 [1 - n L {T NN )} 2 . 



(11) 
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Since the relative phases of the different contributions are not known we cal- 
culate the cross section of the reaction pn — > dM as the incoherent sum 

<T pn -+ dM = N[aW+aW + „(.«)] . (12) 

The dash-dotted curve in Fig. 2 represents the sum + + without 
correcting for the ISI. This result clearly overestimates the empirical cross 
section. In order to achieve agreement with the data for the reaction pn — > dn 
at Q > 10 MeV the theoretical cross section has to be multiplied by a nor- 
malization constant N of 0.68 (see the solid curves in Fig. 2). It is remarkable 
that this reduction factor is not very different from the value obtained from 
the estimation of ISI effect based on Eq. (11). We see that as a confirmation 
of the physics contained in our model being basically correct. 

Comparing the dash-dotted with the dashed curve makes clear that we get 
significant contributions from the exchanges of vector mesons. We want to 
emphasize that those contributions are primarily due to the p exchange. The ou- 
exchange term is negligibly small and the result without its contribution would 
be indistinguishable from the full result shown in the figure. As a matter of 
fact, the magnitude of the p-exchange term is even larger than the 7r-exchange 
contribution. This is not unreasonable. Actually, recently C. Wilkin has argued 
[7] that a large p exchange might be favourable for explaining the observed 
large cross section ratio <y pn ^ pm J a pp ^ ppri [19]. We would like to stress that the 
dominance of the p exchange over n exchange in our model is not due to the 
soft form factor used for the latter. Even with a value of A„- = 1.3 GeV as 
used in the full Bonn model p exchange would still dominate. 

The predictions for the reaction pn — > drj are shown in Fig. 3. The dashed 
curve shows the cross section when only n exchange is taken into account. The 
solid and dash-dotted lines include the contributions from p and u exchange. 
Like in case of n production the u exchange gives only a rather small contri- 
bution. Thus the difference between the dash-dotted and the dashed curves is 
entirely due to the p-exchange contribution. Since there is still a large uncer- 
tainty in the experimental value of the 7p — > n'p cross section (1 pb - 1.5 pb) 
[14] on which the p exchange contribution is based we present calculations 
employing the lower and upper bounds. The corresponding results (lower and 
upper solid and dash-dotted curves in Fig. 3) reflect this uncertainty in the 
input of our model. In any case, the p exchange increases the cross section of 
the reaction pn — > dn' by almost a factor of 20-40. This strong enhancement 
is caused by the rather large pN — > n'N S-wave amplitude extracted from 
the experimental 7p — > n'p cross section. There are indications that the latter 
reaction is dominated by an S , n(1890) resonance, cf. Ref. [14]. 

Admittedly, the ambiguity involved in the p exchange introduce fairly large 
variations in the predicted cross section. But still one can conclude from our 
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model calculation that the cross section for the reaction pn — > n'd for Q ps 30 
MeV should be about 1 ub or even more, and therefore should be measurable 
at new accelerator facilities such as COSY in Jiilich. 

Note that the results shown by the solid curves in Fig. 3 were obtained by 
using the same reduction factor 0.68 for the ISI as in the reaction pn — * dn. 
Estimate (11) suggests a slightly smaller reduction of the pn —> dn' cross 
section (Aisi = 0.83). 

In Fig. 3 we show also the data for the reaction pp ppn' which are from [20] 
(open circles) and [21] (filled circles) for comparison. Evidently, near threshold 
the cross section predicted for the reaction pn — > dn' is significantly larger 
than the cross section of the reaction pp — > ppn' . This is similar to the case 
of n production near threshold (see, e.g. [1]) and can be explained by a more 
favourable isospin factor as well as by the larger phase space in the reaction 
pn — > dM as compared with the reaction pp — > ppM. 

It can be seen from Fig. 2 that the experimental cross section for very small 
excess energies lies above the calculated cross section. In order to demonstrate 
this more clearly we have divided the data points by the result of our model 
(solid curve in Fig. 2) and present this ratio in Fig. 4. Then the enhancement 
of the cross section close to threshold is rather obvious. It is generally believed 
that this enhancement of the cross section at small excess energies is due to a 
strong and attractive interaction of the n meson with the nucleons in the final 
state [2] (a similar enhancement has been also found in the reaction pp — > ppn 
[22]). Therefore, in the following we want to discuss the effect of this final-state 
interaction (FSI) on the reaction pn — » dn. 

Evidently the enhancement occurs only for Q < 10 MeV which corresponds 
to characteristic momenta of p 3 = p^ m below 50 MeV/c. This indicates that 
the enhancement should be mainly due to the coherent final-state interaction. 
Therefore one can safely assume that the range of the FSI is much larger than 
the range of the interaction that is responsible for the production of the n 
meson. In this case the (S-wave) production amplitude and the (S-wave) FSI 
can be factorised [23] and the FSI can be taken into account by multiplying 
the production cross section calculated above with a so-called enhancement 
factor [23] resulting from the long range rjd interaction. 

In the simplest version of this approach the enhancement factor depends only 
on the 77c? scattering length A vd (see, e.g., Ref. [24]). Values for this scat- 
tering length can be readily found in the literature [25,26]. In the present 
investigation, however, we adopt the original prescription of Goldberger and 
Watson [23] where also the nd wave function in the continuum, ip P3 (r), is re- 
quired. Furthermore, we apply two different approaches for the calculation of 
the enhancement factor for the reaction pn — > dn: The optical potential (OP) 
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method and the Foldy-Brueckner adiabatic approach based on the multiple 
scattering (MS) formalism (see also Ref. [23]). Note that both these methods 
have already been used for the calculation of the enhancement factor for the 
reaction pd — > 3 He rj (see Refs. [24] and [27], respectively). 

Let us first describe the OP method. The optical potential in momentum space 
can be expressed through the single ^-scattering term 

V opt (q) = -— — ^ — [tnp(k VP ) + t vn (k vn )]p d (q) ■ (13) 



Here \x is the rjd reduced mass, and t V N the r/N t-matrix which is related to 
the scattering amplitude f V N by t^ik^) = (1 + ^) f V N{k n N)- Note that 
we use the scattering length approximation for the latter, i.e. f V N(k V N) = 
((ctriN)' 1 — ik V N) 1 - krjN is the modulus of the relative r/N momentum and 
Pd(q) the deuteron form factor, 

Pd(q) = J d 3 r exp ^|f) \V d (r)\ 2 . (14) 

— * 

In the r]d cm. system k V N is given by 

r m m v Pd m + m v /2 ^ 

k v n= : P v : = : Pa- (15) 

m + m v m + m v 2 m + 



The optical potential (13) is transformed into coordinate space and then the 
Schrodinger equation is solved numerically yielding the continuum wave func- 
tion ippzif^d) as well as the ^-scattering length. The enhancement factor is 
simply given by \%(p 3 ) = |^ P3 (0)| 2 [23]. 

In the multiple scattering formalism we use the Foldy-Brueckner adiabatic 
approximation and start from the rjd wave function defined at fixed coordinates 
of the proton (r p ) and the neutron (r n ) (see Ref. [23] for details): 



ipk(r v , r p , f n ) = exp (ikr v ) 



where 



+ £ (exp («>„) + t m e -*^TA exp { f kf n) ) 

t vn exp(ikr n ) ( - exp (ikr pn ) ,.f^\ 
+ -=r exp (ikr n ) + t vp exp (ikr p ) , 



D=U-t^ e ^^). (16) 

V pn J 
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Here f pn — f p — r n and k = p 3 and k, r pn , etc., are the moduli of these vectors. 
t V N is defined in Eq. (13). Then the 7/ci-scattering length and the enhancement 
factor follow from 



^ms _ m d / 2t ?? jy(A: r? jv — 0) 
vd m v + m d 1 - t vN {k vN = 0)/r 
A^Kps) = I (lUr, = 0, f p = r72, f n = -f/2)> | 2 . 



The resulting enhancement factors Afsi(P3) are plotted in Fig. 4 (upper panel 
for the multiple scattering formalism and lower panel for the optical potential) 
for different values of the T/iV-scattering lengths given in Table 1. In this Table 
we also present the calculated values of A®J and and compare them with 
the results of Ref. [26]. Of course, the optical potential model is a questionable 
approximation for the r/d system. It only gives reasonable values for the rjd 
scattering length if a V N is small so that the first iteration is dominant. For 
larger a V N it might still give a reasonable estimate of provided that there 
is a cancelation between different higher-order iterations. 

Nevertheless, it is expected that the multiple scattering approach gives a more 
reliable estimation for A v d- It is known that the Foldy-Brueckner approach 
describes the 7r<i-scattering amplitude at low energies quite well (see, e.g., 
Ref. [32]). The rj- meson is certainly somewhat heavier than the pion but its 
interaction with the nucleons is likewise considerably weaker than the forces 
in the NN system. Therefore the adiabatic approximation might still be quite 
reasonable also in the case of r/d scattering. For example, we see from Ta- 
ble 1 that within the MS approach an attractive r]N interaction never yields 
spurious repulsive results for the r]d system — as it happens for the optical 
potential and for the first approach used by Green et al. [26]. We want to point 
out that our results are in rather good agreement with the second approach 
of Ref. [26] for a wide range of a V N values. 

The solid curves in Fig. 4 show our results based on the 7]N scattering a V N = 
0.291 +i 0.360 fm derived from r\ photo-production data [29]. Evidently, they 
describe the enhancement of the cross section rather well. The dashed and 
dotted curves correspond to the values a V N = 0.3 + i 0.3 fm and a vN = 0.476 + 
% 0.279 fm proposed in Refs. [24] and [27], respectively. The former curve which 
corresponds to the smaller value of a V N still agrees satisfactorily with the data, 
while the latter one overshoots them. The dash-dotted curve based on a vN = 
0.717 + % 0.263 fm from a recent analysis of Batinic et al. [31] overestimates 
the data at Q < 10 MeV by more than a factor 2. Thus, we conclude that our 
analysis of the reaction pn — > dr] close to threshold favours fairly small r)N- 
scattering lengths, Re a V N < 0.3 fm. The t]d scattering length corresponding 
to these values are A rj d ~ 0.4 + i 1.2 fm, cf. Table 1. We would like to point 
out that the above r/N scattering lengths are considerably below the lowest 
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value, Re a V N ~ 0.7 fm, for which an rjd bound state is expected to exist [26]. 
Finally, we want to call attention to the fact that both employed approaches 
give very similar results for the enhancement factor. Furthermore, we would 
like to remark that a recent Faddeev-type calculation of the rjd system testifies 
that the MS approach is quite reliable for small values of Re a V N like those 
supported by our analysis [33]. 

In summary, using the two-step model we calculated the cross sections of the 
reactions pn — > drj and pn — > drj' near the threshold including n, p as well as 
lo exchange. The u contribution is found to be rather small. The contribution 
of the p exchange is of comparable magnitude to the one from 7r exchange in 
the case of rj production and might be even more important for rj' production 
due to the presence of a possible 5 , n(1980) resonance for which the SAPHIR 
collaboration found strong evidence in the photo production of n' mesons near 
threshold. Our results agree with the experimental data on n production if 
we take into account effects from the initial state interaction. We considered 
also the influence of the final-state interaction using the optical-potential and 
multiple-scattering approaches. Both methods yield similar results for the en- 
hancement factor. Agreement with the data could be obtained only for rjN 
interactions with small ^A-scattering lengths, Re a V N < 0.3 fm, i.e. values 
that are considerable below the lowest value, Re a V N ~ 0.7 fm, for which an 
rjd bound state is expected to exist. 

We are grateful to H. Stroher for useful discussions. 
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Table 1 



77-deuteron scattering lengths in fm resulting from the optical potential method 
(OP) and the Foldy-Brueckner approach (MS) as described in the text. 'Greenl' 
and 'Green2' refer to results presented in Ref. [26]. 







A MS 

rjd 


^Grccnl 
rjd 


^Grccn2 
r)d 


Bennhold-Tanabe [28] 


0.25 + 


i0.16 


0.606 + i0.89 


0.613 + iO. 602 


0.66 + iO. 71 


0.66 + i0.58 


Krusche [29] 


0.291 - 


- i0.360 


-0.14 + il. 33 


0.399 + il. 18 


0.17 + il. 35 


0.42 + il. 25 


Wilkin [24] 


0.30 + 


i0.30 


0.08 + il.37 


0.526 + il. 06 


0.39 + il. 28 


0.58 + il. 11 


Faldt [27] 


0.476- 


- i0.279 


-0.15 + i2. 32 


1.01 + il. 387 


0.81 + i2.15 


1.22 + il. 56 


Batinic [30] 


0.888- 


- i0.274 


-2.06 + il. 37 


1.85 + i2.83 


-2.90 + i4.12 


2.37 + i5. 79 


Batinic [31] 


0.717- 


- i0.263 


-1.94 + i2. 31 


1.62 + i2.11 
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a) b) 
^\>v/\/^r|(r|') p — « "Y^>^ m(m') 

71 * x , 71 * V 



d 

P 



c) d) 



V\>v/\/^ t|(T|') n — » v \>v/v^ T|(ri') 

Ji°/\n ity\p 
— ^ W= d p , / \ a d 



Fig. 1. Diagrams describing the two-step model (TSM). Note that besides the 
7r-exchange contribution also diagrams involving the exchange of the p and to mesons 
are taken into account. 
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A/s-m 3 -m 4 (MeV) 



Fig. 2. Cross section of the reaction pn — > drj as a function of the cm. excess energy. 
The dashed curve shows the result from the 7r-exchange contribution alone whereas 
the dash-dotted curve is the sum of ir, p, and u> exchange. The solid curve represents 
the results including all contributions (it, p, oS) multiplied with a normalization 
factor N = 0.68 in order to take into account effects from the initial state interaction 
(see text). The data points for are taken from Refs. [1] (open circles) and [2] (filled 
circles) . 
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10 20 30 40 50 60 70 

A/s-m 3 -m 4 (MeV) 

Fig. 3. Cross section of the reaction pn — > dr/' as a function of the cm. excess 
energy. The dashed curve shows the result from the 7r-exchange contribution alone 
whereas the dash-dotted curves are the sum of it, p, and lo exchange. The solid 
curves represent the results including all contributions (ir, p, u) multiplied with a 
normalization factor N = 0.68 in order to take into account effects from the initial 
state interaction (see text). The upper and lower solid and dash-dotted curves are 
the results obtained using the maximal and minimal values of the p°p — ► rj'p S-wave 
amplitudes. The data points are from Refs. [20] (open circles) and [21] (filled circles), 
respectively. 
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pn^dr| 

A. FSI (multiple scattering) 



T 



pn— >dr| 

A. FSI (optical potential) 



" " " * ~ ~ ~~- -*-- r - -~ :T ' ~ * „ 



3 
2 
1 


5 10 15 20 25 30 

Vs-m 3 -m 4 (MeV) 

Fig. 4. The enhancement factors Afsi(P3) for different values of the r/A^-scattering 
lengths as a function of the cm. excess energy. The solid curves corresponds to the 
r]N scattering length a V N = 0.291 + i 0.360 fm derived from 77 photo-production 
data [29]. The dashed and dotted curves correspond to the values a v N = 0.3 + z0.3 
fm and ci v n = 0.476 + i 0.279 fm proposed in Refs. [24] and [27], respectively. The 
dash-dotted curve is the result for a V N = 0.717 + i 0.263 fm, found in a recent 
analysis of Batinic et al. [31]. 
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